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Q ■ Abstract 

We explore the possibility that the observed baryon asymmetry of the universe 
is the result of an earlier phase transition in which an extended gauge sector breaks 
down into the SU{3) C x SU(2) L x U(1) Y of the Standard Model. Our proto- 
typical example is the Topflavor model, in which there is a separate SU{2)\ for 
the third generation from the SU (2)2 felt by the first two generations. We show 
that the breakdown of SU{2)\ x SU (2)2 — ► SU{2)l results in lepton number being 
asymmetrically distributed through-out the three families, and provided the SM 
electroweak phase transition is not strongly first order, results in a non-zero baryon 
number, which for parameter choices that can be explored at the LHC, may explain 
the observed baryon asymmetry. 
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1 Introduction 



The origin of the observed baryon asymmetry of the universe is one of the most fundamen- 
tal open questions in particle physics and cosmology. Numerically, the baryon-to-entropy 
ratio is 

Vb = y = 9 - 2 ~o" x 10_n > (!) 

where % = — ni and s are the baryon number density and entropy, respectively. 
Famously, Sakharov has shown [1] that in order to generate a baryon asymmetry in an 
initially baryon-symmetric universe, there must be: (1) baryon number (B) violation; (2) 
C and CP violation; and (3) a departure from thermal equilibrium. These requirements, 
especially in light of null results obtained by experimental searches for S-violation (for 
example, proton decay), are difficult to realize in the framework of an electroweak (EW) 
scale model, prompting attention to high scale mechanisms such as GUT baryogenesis 
[21 O HI [5] , leptogenesis [SIE] and the Affleck-Dine scenario [8]. Such mechanisms can be 
viable, but are difficult to test experimentally. 

However, the Standard Model (SM) already contains a way to reconcile large baryon- 
number violation at the EW scale with the lack of experimental evidence for such inter- 
actions at low scales. As noted by Kuzmin, Rubakov and Shaposhnikov [9], the baryon- 
number violation induced by sphaleron transitions between two different iV-vacua of the 
SU (2) l electroweak interaction are unsuppressed at temperatures larger than the EW 
scale, but are so highly suppressed at low energies as to be essentially irrelevant. This 
idea of electroweak baryogenesis (EWBG) [10J is very attractive, and could have been re- 
alized within the SM. EWBG makes very specific demands of the electroweak symmetry- 
breaking (EWSB) sector of the theory, and leads to observable consequences at future 
colliders. For example, the need for S-violating interactions to be out of equilibrium at 
the time of the phase transition requires that the electroweak phase transition is strongly 
first order, with 

^ > 1. (2) 

J- c 

This in turn puts constraints on the Higgs potential, and demands a small Higgs quartic 
(and hence Higgs mass). 

Electroweak baryogenesis in the SM is thus ruled out; the LEP-II bound on the 
SM Higes mass [11] , m^ M > 115 GeV is incompatible with a the requirement that the 
sphaleron processes be out of equilibrium at the phase transition, Eq ([2]). While this ren- 
ders the SM itself unable to produce the asymmetry, it illustrates the fact that experiments 
at the TeV scale are able to directly probe EWBG. Extensions of the SM which include 
new physics at the EW sea can in fact re-open the possibility of EW baryogenesis by 
introducing additional bosonic fields [TS1 EH HH US, HE] as in the Minimal Supersymmetric 

1 In addition to new fields, a greatly increased Hubble expansion[27l [28l [29l [30] can re-open some of 
the Higgs mass parameter space. Such models usually require the very fast Hubble expansion to only 
happen close to the electroweak scale. 
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Standard Model (MSSM) [HllIEllISll2nil2Ill22ll231l23or fermions strongly coupled to 
the Higgs [311 E2]. Eventually, future colliders such as the Large Hardron Collider (LHC) 
and possibly International Linear Collider (ILC) will unravel the nature of EWSB, and 
should provide a better understanding as to the nature of the electroweak phase transi- 
tion (EWPT). If it proves first order, this will be a crucial piece of indirect evidence for 
EWBG. If it does not, it will raise the interesting question: Does the idea of electroweak 
baryogenesis have to be abandoned? 

While the lack of a first order EWPT would strongly disfavor electroweak baryogenesis, 
it may be that the basic picture of baryon-number violating processes arising from non- 
perturbative gauge theory dynamics, impotent at low energies but unsuppressed at the 
TeV scale, can survive. Since many extensions of the SM predict the existence of new 
non-Abelian gauge symmetries at the TeV scale, it may be that some theories already 
addressing unrelated problems may in fact contain the ingredients necessary for "electro- 
weak style" baryogenesis. A new gauge group G could break down somewhere above the 
EWPT, and generate an asymmetry through its own strongly first order phase transition. 

This is a novel idea, but one which is not entirely straight-forward to realize in a 
realistic setting. Below this new phase transition, the theory is still SU(2)l X U(1)y 
symmetric, and the usual EW sphalerons are unsuppressed. They will try to wash-out 
any generated baryon asymmetry unless B — L ^ 0. In principle, one could arrange the 
representations of the SM matter under the new gauge symmetry G such that B — L is 
not preserved due to the nonzero mixed anomaly G-G-U(1)b-l, but this kind of non- 
universal assignment of representations is also somewhat delicate. The requirement that 
G-G-U(1)y gauge anomalies cancel will in general require more chiral fermions charged 
under G and the SM. Thus, we will not consider this possibility in detail. Instead, we 
employ a more subtle mechanism in which G acts only on the third generation fermions, 
and thus B — L is conserved (within a generation). The breakdown of G can thus be 
accompanied by production of an asymmetry among the third family quarks and leptons. 
Baryon number will quickly equilibriate among the three families because of the large 
quark masses and mixings, but each lepton family number will remain distinct because 
of the tiny neutrino masses. Above the EWPT, the EW sphalerons will result in B = 
and L\ = L 2 = — L 3 /2 ^ 0. Provided the EW sphalerons are in equilibrium during the 
EWPT, as the fermion masses turn on from EW symmetry-breaking, a non-zero (though 
diluted) baryon number will result [3"3"l IM] . So in fact this mechanism requires that the 
EWPT not be strongly first order. 

The outline of this paper is as follows. In Section [21 we review how a baryon asymme- 
try can be generated even when B — L = 0, provided there is a non-zero asymmetry in the 
third family lepton number, and how this eventually translates into a baryon asymmetry 
after the EWPT. We apply this mechanism to the "Topflavor" model [35] which is known 
to contain non-perturbative interactions which violate baryon- and lepton-number in the 
third family [36J. In fact, we find that it is possible to generate a baryon asymmetry of the 
right magnitude. The Topflavor model, phase transition, CP violating sources, diffusion 
equations and calculated baryon number density are considered in Section [3J Our results 
show that the right baryon number density can be generated for parameters that would 
render this model testable at the LHC. In Section [3] we present our conclusions. 
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2 B from a Family Asymmetric Distribution of L 



In this section we review the mechanism by which an initial condition that has B = L = 
can nonetheless result in a non-zero baryon number density provided L\ = L 2 = — L 3 /2 7^ 
[551 S3- We will show how the specific example of the Topflavor model can generate 
these initial conditions through the non-perturbative dynamics of its phase transition (in 
a very similar way to that in which baryon number is generated in a traditional EWBG 
scenario) in Section [3j Thus, for now we assume that the initial phase transition has 
generated L3 = B = A and L\ = L2 = 0. The unsuppressed EW sphalerons will rapidly 
evolve this into a state with B = L = but L\ = L2 = —A/3 and L3 = +2A/3. We 
now study what happens as the universe moves through the EWPT, and show that the 
non-zero r lepton mass will result in B ~ 10~ 6 A. Our discussion closely follows that of 
Ref. [33]. 

The smallness of observed neutrino masses indicates that lepton-number violation is 
out of equilibrium at the EW scale. Thus, each lepton flavor has a seperate chemical 
potential fii with i = 1,2,3. This is in contrast to baryons, because the large quark 
masses and mixings keep baryon flavor violation in equilibrium, and thus the quarks are 
described by a single chemical potential fi. We consider the SM matter consisting of three 
families each of which consists of two quarks (an up- type and down- type) with masses m qi , 
a charged lepton of mass m;. and a neutrino which for our purposes can be approximated 
as massless. The free energy per unit volume for the system in equilibrium at temperature 
T is given by 

6 3 

T = 6 F K./^) + YP F ( m i*^0 + ' ( 3 ) 

i=l i=\ 

where the SU(2)l gauge interactions maintain equilibrium between the charged lepton 
and its neutrino and the up- and down-type quarks of a given generation. The free energy 
density for a (single helicity of a) fermion of mass m and chemical potential /x is given by 

F(m,fi) = -T J ^[ln(l + e-^/ T )+ln(l + e-^/ T )], (4) 

where Ki is the spatial momentum of the fermion, and E = Kf + m? is its energy. At 
high temperatures, T ^> m, /i, this may be approximated as, 

F (m, /,) » F (m, 0) - ^ 2 T 2 (l - JL^ . (5) 

The (individual) leptonic and baryonic number densities may written 

U = —[2F(m,fH)+F(Q,iH)]n~tiiT*pi, (6) 
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where 

3 ml. 1 mf 



a: 



2vr 2 T 2 ' ^ vr 2 T 2 

i=l 



Electroweak Sphaleron transitions violate ^2 Li and -B, but preserve the three combi- 
nations Aj = Li — B/3. In terms of the chemical potentials these are 

A^Li-iBw^a-^A, (9) 

We can invert the above relations to obtain each /ij in terms of the quark chemical potential 
fi, temperature T, and the conserved value of the corresponding Aj. Effectively, the EW 
sphalerons convert nine quarks and one lepton of each family into nothing. In thermal 
equilibrium, this leads to the relation \x = — Yli^i/^- Using this fact, together with the 
three conservation equations Eq (|2J), allows us to express the quark chemical potential in 
terms of the values of the Aj, 

which can be combined with Eq (JBJ) to obtain the final baryon number density [33J 

b=i ir } r (11) 



-WE^V B-L = 0. 

The first of these results is the familiar relationship applicable to theories that directly 
generate a non-zero B — L (such as leptogenesis) and indicates that in such theories 
primordial B cannot be completely washed out, and a primordial L will be converted into 
B by EW sphalerons. The second result shows how in a theory with B = L = but 
the individual Aj non-zero, the turn on of the charged lepton masses will also generate a 
non-zero B. In the scenario we are considering, with initially B = and L3 = 2 A/3, and 
taking the freeze-out temperature to be the close to the EW scale, the resulting baryon 
number is diluted to about B ~ 10~ 6 A [331 E3j- 

Since the dilution factor plays a relevant role in our work, let us expand on its origin: 
To compute the above quoted dilution factor, we have assumed a second order phase 
electroweak phase transition. Under this condition, the sphaleron processes will remain 
in equilibrium until the weak spahleron rate is of the order of the expansion rate of the 
Universe. The departure from equilibrium therefore occurs at the freeze-out temperature 
Tp, such that v(Tf)/Tf — 1. Using the relation m T (T) ~ h T /y/2v(T), and the condition 
v(Tf)/Tf = 1, we get that the final baryon number is approximately given by 

4 h 2 

B~ -A-^ ~ -1.6 x l(r 6 A. (12) 

13tt 2 2 V ; 
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3 Lepton Number Generation in the Topflavor Model 



3.1 The Topflavor Model 

The gauge extension of the SM that we consider is based on the gauge group SU(3) C x 
SU{2)i x S77(2) 2 x U(1)y- While the SU(3) C and U(1)y subgroups remain the same as 
those of the SM, the SU(2) L group of SM is expanded to a larger SU(2) 1 x SU(2) 2 in a 
flavor dependent way. The fermion content of the model is identical to SM. Under the new 
SU (2)i and SU(2)2 groups, the doublets of the third generation transform as doublets 
under 377(2) i and singlets under SU{2) 2 , while the first and second generation doublets 
are singlets of SU(2) 1 and doublets of SU{2) 2 . Thus, their SU(2) 1 x SU{2) 2 x U{l) Y 
quantum numbers are 

Q 3 = (2,l) 1/6 , Q 1 ' 2 = (l,2) 1/6 , 

L 3 = (2,l)_ 1/2 , L 1 * = (l,2)- 1/2 . (13) 

After symmetry breaking, the Standard Model SU (2) l group emerges as the unbroken 
diagonal subgroup of SU '(2) 1 x SU(2) 2 . The corresponding SU(2) L gauge coupling is 

9L = (14) 
V9i + 92 

This relation implies that when one of the gauge couplings becomes large, the other one 
approaches gi from above, and thus both g\ and g 2 are necessarily larger than g^. Thus, 
a convenient parameterization is given by, 

9i = ——, 92 = T, (15) 

sin (p cos <p 

in terms of an angle 0. We will use a simplified notation s = sin and c = cos <f> below. 

The symmetry breaking of the extended gauge group, SU (2)i x SU (2) 2 to the SM weak 
gauge group SU (2) l is accomplished by introducing a vacuum expectation value (VEV) 
to a scalar field S, which transforms as a bidoublet (2, 2) under the extended gauge group 
transformations. After the S77(2)i x SU(2) 2 breaking, E can be decomposed under the 
residual diagonal SU(2)l symmetry into a complex singlet a and a complex triplet r (half 
of which is eaten by the Z' and W's) , 

E=if a + T3 ^ I . (16) 



2 

We introduce the scalar potential, 



V2- 



r_ a - r 3 



Ve = m 2 |S| 2 + A|(SS)| 2 + A / |S| 4 +(-i/J(SS) + A(SS)|S| 2 + / i .c. ), (17) 
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where D and A are complex parameters and we use a notation that supresses the gauge 
indices: (EE) = E a 5E c je ac e5j and |E| 2 = Tr(EtE). For appropriate choices of parameters, 
this potential results in the VEV, 

(Stf> = lu e i6 °5 fk . (18) 

which will generally be complex, and will be solution of the equations, 

2 De 2i0 ° + D*e- 2l0 ° - m 2 

U ° ~ A + A' + Ae 2 ^ + \* e - 2W ° 



— acos Re 
4 



-2D* + X*u 2 
-2D + \u 2 



(20) 



Choosing some representative parameters, taking m = 200 GeV, D = 5 x 10 5 e* GeV 2 , 
and A = A' = A = 0.05, we obtain a zero temperature VEV described by uo — 2.7 
TeV and 6q — —0.7. This particular set of parameters has been chosen with small quartic 
interactions in order to have a first order phase transition, with the dimensionful quantities 
arranged such that the SU (2) x SU (2) breaking scale is of order TeV. Precision electroweak 
constraints have been extensively considered in the literature [35], and typically require 
u > a few TeV. Requiring that the extended instantons of the strongly coupled 577(2) 
do not mediate unacceptably large proton decay [36] further requires the gauge couplings 
to satisfy sin 2 (p > 0.2. We will illustrate our results with the representative point chosen 
above, and s 2 = 0.4, in order to be consistent with all constraints. 

The fermion doublets of either SU(2) 1 or SU (2) 2 transform as doublets under SU{2) L . 
In addition, there is a SU (2) L triplet of heavy gauge bosons from the breaking. We denote 
the neutral and charged heavy gauge bosons as Z' and W r± . Their masses are degenerate 
and given by 

M w , ± = M z , = {gl + g\)u 2 = ^u 2 . (21) 

At the electroweak scale, v ~ 174 GeV, the remaining SU(2)l x U(l)y electroweak 
symmetry is broken to U(l) em as in the SM. This is accomplished by giving a VEV to one 
or more Higgs boson doublets. There are two possible representations for the Higgs boson 
under SU{2)i or 5£/(2) 2 , either (2, l) ±1 / 2 or (l,2) ±1 / 2 - We will focus on the first case 
(sometimes called the heavy case) as it motivates the large third family fermion masses as 
they are the only family whose Yukawa interactions are 5?7(2)i x 5f/(2) 2 gauge-invariant. 
In a non-supersymmetric theory, a single Y = +1/2 Higgs doublet H u suffices, but our 
results are largely unchanged if we generalize to a Y = ±1/2 pair of doublets, H u , H d 
instead. In order to connect more easily with the supersymmetric case, we consider the 
case with two Higgs doublets below. 

The Yukawa couplings for the first two generations can be generated by adding an 
additional "spectator" Higgs-like doublet H' u (in a supersymmetric theory it would be 
a pair of doublets including H' d ) charged under 5?7(2) 2 . They couple to the first two 
generations via Yukawa couplings and mix (slightly) with the regular Higgs(es) via inter- 
actions such as A 1 H u T l H' d . The small Yukawa couplings for the first two generations can 
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be naturally obtained from the small mixing angle between H' and H . The full set of this 
kind of gauge- invariant H-H-H' interactions between Higgs(es), spectator Higgs(es), and 
£ include, 

A x HJlH' d + h.c; A 2 H d Y J H' u + h.c, 
A[H u ^H d + h.c; A' 2 H d Y)H' u + h.c, 
c^H'JlHl + h.c; c^'H' u T)Hl + h.c, 
Cl n*H d Y,Hl + h.c; c^n'H'jpHl + h.c (22) 

These interactions will be important below, because they (indirectly) drive the Topflavor 
phase transition's production of L 3 . 



3.2 The Phase Transition 

The details of the SU {2)i x SU{2) 2 — > SU(2) L phase transition depend on the finite 
temperature effective potential. In order to generate a baryon asymmetry in the Topflavor 
model, the processes which violate baryon- and lepton-number must be out of equilibrium 
at the time of the phase transition. We are interested in the regime where g\ ^> g 2 , for 
which we can approximate the sphalerons associated with the SU(2)i x SU{2) 2 symmetry 
breaking as purely arising from SU (2)i. Their rate, 

u 7 

T sph ~ k — Exp [—E sp h/T] , (23) 

(with E sp h = inu/gi and k a dimensionless parameter of order one |37J) must be much 

less than the Hubble expansion, H ~ g X J 2 T 2 /mpi. Assuming k ~ 1, g* ~ 100 and 
u c^T c ~ TeV, this requires, 

" > 2.5, (24) 



9iT t 

and for s 2 ~ 0.4, we should hav^l u c /T c > 2.5. 

The potential for S, Eq. f fTTl) . can be expanded to 

1 




V (a) = |m>| 2 + — —\n~\ '. r . j —^D(t'~ + /'•'•• 1 (25) 

where for brevity we have not shown the terms involving the triplet component r. At one 
loop, there are both temperature-dependent and temperature-independent corrections to 
the potential, 

V(u,9,T) = V (u,9,0)+Vi(u,9,0)+Vi(u,9,T) . (26) 



2 Note that such a strong first order phase transition may provide a strong signature in gravitational 
waves [40l IHl |42j H31 HI] detectable at the planned space interferometer, LISA. We will pursue this idea 
in a separate paper [45]. 
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We consider the limit where the gauge couplings of the SU (2)'s are much larger than 
that of the E self-interactions A and thus approximate the complete one-loop corrections 
by those from the gauge bosons, Z' and W. We renormalize parameters according to a 
scheme that absorbs the correction to the position of the zero-temperature minimum, and 
find, 
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647T 2 \ S 2 C 2 



2 \ 2 
9t \ 2 

L 1 u 



and finite temperatue correction, 



u 



lo s — 

o 



u 



2 J + 2m o 



(27) 



2tt 2 



T 



with 



(p!f±) = J™ dx ■ x 2 {log (l - e -V*W/( S W)) _ lQg(1 _ e - x) j 



(28) 



(29) 



where $ = 9. The Debye screening effect on the longitudinal modes of massive gauge 
bosons is neglected because their corrections g{F is small compared to their masses in- 
duced by the Higgs vev[47J. Note that the one loop corrections from the gauge sector 
depend only on the magnitude of the VEV u and not on its phase. 



For the sample parameters chosen above, 



rn 



200 GeV, D = 5 x 10 5 e l GeV^ 



A = A' = A = 0.05, and s 2 = 0.4, we find that the critical temperature for these parameters 
is T c ~ 840 GeV, and the VEV at T c is described by u c ~ 2.7 TeV, 9 C ~ -0.7, indicating a 
first order phase transition that is easily strong enough. In Figure [1] we plot the effective 
potential for several choices of temperature. 



3.2.1 Bubble Profile and Evolution 

At T c , the rate for nucleation of bubbles with (a) ^ becomes large, and the bubbles 
expand to fill the universe with the true vaccuum. In this subsection we make some rough 
estimates of the properties of the nucleated bubbles, which are pertinent to the eventual 
generation of baryon asymmetry as they provide the out-of-equilibrium dynamics which 
results in lepton-number being unequally distributed through-out the three generations. 
We will simplify the treatment by considering the phase transition as a quasi-equilibrium 
process such that the temperature is varying slowly enough that the properties of the 
nucleating bubbles can be obtained by studying fixed temperature solutions at T ~ T c . 

Assuming that the nucleated bubbles have spherical symmetry [38], the Euclidean 
action of the configuration become^! 

S 3 (T) = 4vr J 'c*rr 2 {(v(a>) 2 + V((a),T c )} (30) 
3 We use the 0(3) approximation since the nucleation temperature T n ~T C > 1/(2R) 46 . 
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u (TeV) 



Figure 1: The effective potential as a function of the magnitude of the VEV for three 
different choices of temperature. The phase of the VEV at each point is chosen as the 
solution of the Equations of motion for that value of the magnitude. 



where (a) is a function of radius r that describes the configuration. The transitions will be 
predominantly mediated by the solutions which minimize this action, give by the solutions 
of the equations, 

d 2 u 
dr 2 

~d 2 6 2d6~ 
+ — r- 



2du 
r dr 



u 

subject to the boundary conditions, 

du 

dr r=0 

dd 



dr 2 



r dr 



6V 

Su 
5V_ 



(31) 



0; 



0; 



u 



0\ 







u=0 



(32) 



These coupled differential equations are somewhat difficult to solve. Instead of looking 
for detailed solutions, we will use an ansatz for the profile, and a variational approach 
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to determine the parameters that describe it. We write the solutions in the form of two 
"kinks" with the proper asymmptotic behavior, 

u{r) = [1 — Tanh (a (r — R))] 

B{r) = 6 U=0 + ° c ~° u=0 [1 - Tanh (a (r - R))} (33) 

where the radius of the bubble is ~ R and the width of the bubble wall is ~ 1/a. In 
addition to imposing the form of the solution, we also assume that the bubble width is 
the same for u and for 9 (which we expect to be approximately true). We determine a 
numerically by plugging the solutions of Eq. (13"3"|) into the action Eq. (|3"Uj) and finding 
the value of a which minimizes the action. For our standard parameter choice we find 
a ~ T/2 and L w = 1/a ~ 2/T. The profiles are plotted in Fig [2j We expect that for 
large bubbles, the details will become independent of R, which in fact proves to be true. 




r (TeV 1 ) 



Figure 2: The bubble profile at the critical temperature T c for R = 10 TeV . 

Under our quasi-equilibrium assumption, the expansion of the bubble is driven by the 
fact that the gauge bosons acquire masses inside the bubble, and thus the free energy is 
minimized for large bubbles [171 [39] . By analogy with the SM EW phase transition, we 
estimate the bubble wall velocity v w ~ 0.05, though we find that our final results are very 
insensitive to the precise value of v w . 
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3.3 Diffusion Equations in the Topflavor Model 



We now compute the prediction for the L 3 generated during the transition in which 
the Topflavor model breaks down to the SM. The underlying picture is similar to the 
standard EWBG picture in the SM (or MSSM). The bubble of true vacuum is expanding, 
and generates chiral charge through the CP- violating interaction of the plasma with the 
bubble wall. In the specific case of Topflavor, the particles which interact strongly with 
the wall are the Higgses, through the interactions in Eq f[2"2"j) . These charges diffuse freely 
in the unbroken phase and are converted into B and L3 by a combination of the Yukawa 
interactions and the unsuppressed sphalerons. As they pass into the broken phase, they 
are frozen. 

In the limit g\ ^> g%, we neglect the SU(2)2 sphalerons associated with the first two 
families. The quark Yukawa interactions and the QCD instantons, together with the fact 
that all of the quarks diffuse at approximately the same rate, allows us to constrain the 
light quark densities in terms of the right-handed bottom density b, 

Qil = Q2L = ~2U R = -2D R = -2S R = -2C R = -2b . (34) 

Thus, the species whose densities we will track are the left-handed top and bottom doublet, 
Q = t L + b L , the right-handed top t = t R , the right-handed bottom b = b R , the left-handed 
lepton doublet L 3 = t l + v T , and the Higgs h = (h+ + h° u — h d — h®). We assume that the 
H-H'-T, interactions Eq. (I2"2"j) are fast enough such that the spectator Higgses H' are kept 
in equilibrium with the Higgs, and thus E = 0, hi = (hf£ + h' u — h' d — h' d ) = h, and we 
do not include the densities S and hi in the diffusion equations. For relativistic particles 
near equilibrium, we can write the number densities in terms of a chemical potential m 
as rii = ki[iiT 2 /6 where hi counts the number of degrees of freedom, 

kg = 6; k L = 2; k t = k b = 3; k h = 8 . (35) 

The diffusion equations will contain the interactions which are fast compared to 
the time scales at which the elements of the plasma are diffusing, n = Di/v w , where 
v w ~ 0.05 is the speed of the bubble wall's expansion and Di is a diffusion coefficient 
which characterizes the interactions with the background plasma. Typically, one expects 
Dq = D t = D b ~ 6/T and D L ~ D h ~ 110/T [49]. Thus, we consider the processes 
characterized by rate r > tq. These interactions include the SU (2)i sphalerons with rate 
r 1; the QCD instantons with rate Tqcd, and the top quark Yukawa coupling to the Higgs 
with rate T y . We continue to assume that the sphalerons associated with S'f/(2) 2 can be 
neglected. These rates are estimated to be equal to [371 150| - 

r„ - yA^ 5 (^yr = 7.4GeV, (36) 

Tqcd ^ Wn'a 4 s T = 0.3 GeV , (37) 
T 1 ~ 30«i 5 T = 0.1 GeV, (38) 

where X t ~ 1 is the top Yukawa coupling, as(T c ) ~ 0.08 is the strong coupling constant, 
and k' ~ 1 is a dimensionless coefficient. We have evaluated the rates for s 2 = 0.4 and 
T c = 840 GeV, as is appropriate for our example parameter set. 
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We approximate the bubble as large, and thus treat the problem one-dimensionally, 
with the z-axis perpendicular to the wall, whose location is at z = 0, with the z > side 
in the broken phase. The rates of change of the various densities are described by the 
coupled set of equation^ 
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v w ti - D h h" 
vj - D Q b" 
v w L'n - DjL'i 



-r,. 



Q_ 

kn 



h t 



6r 



QCD 



Q t b 
2— 9 

kn k t k b 



-6^ 

-r, 
-r 

^QCD 

-2r\ 



3 Q + U 

k Q k L 

kn 



kc 



h 




t 




+ 




k h 


h 


h 




t 




+ 




ku 


h 



+ 3r 



QCD 



Q_ 

kg 



t b 
^-9 — 



+ lh 



Q t b 
2— 9— 

k Q k t k b 

Q Li 

kn k L 



(39) 



where primes denote derivatives with respect to z and 7^ is the CP-violating source for 
the Higgs induced by the bubble wall, approximated as a step function, 



lh 



1h u 1h d 





(-L w <z<0) 
[z > or z < — L,„) 



(40) 



whose magnitude we estimate below. 



3.4 CP-violating Sources from Spontaneous CP violation 



We consider the CP-violation arising from the spontaneous CP-violation associated with 
the phase of the VEV (cr). This field couples directly to the EWSB Higgses H u and Hd 
and thus influences their number densities as they scatter off of the bubble wall. As we 
saw in section 13.21 the phase of (a) varies as one moves from inside the bubble of true 
vacuum to the unbroken phase (see i.e. Figure [2]), 



A9 



6 r . — 6 



u=0 



-acos 



-2D* + \*u 2 c 
-2D + \ul 



+ 



(41) 



and thus is space-time dependent as the bubble expands. 



4 Note that leptons diffuse faster than quarks, and thus the B — L charge density (Q + t + 6)/3 — L is 
not zero locally. 
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Figure 3: Feynman diagrams for the leading contributions to the self-energies of H u in 
the background of a spacetime varying (£). 



In computing the value of the CP-violating source jh, we follow the treatment first 
introduced by Riotto [51], [52] based on the closed time path (CTP) formalism, which 
allows us to capture the main non-equilibrium quantum effects. The CTP formalism 
distinguishes fields with arguments on the positive and negative branches of the closed 
time path. This doubling of fields leads to six different real-time propagators; for a generic 
scalar field <p, 

G>(x,y) = -i{<j>{x)^{y)), 

G<(x,y) = -i(tf(x)<j>(y)), 

G^(x,y) = 6(x,y)G>(x,y) + 6(y,x)G<(x,y), 

Gl(x,y) = 9(x,y)G<(x,y) + 9(y,x)G>(x,y), (42) 

which are conveniently written matrix G(x, y): 

I G t (x,y) -G<(x,y) \ 
G(x,y)= 1 ,U> _ y ,U) . (43) 

\G>(x,y) -G\x,y) J 
From the Schwinger- Dyson equations of the path-ordered two-point functions, one obtains 

Sz / dz [Z>(X,z)G<(z,X)-G>(X,z)X<(z,X) 

J — oo 

+G<(X,z)Z>(z,X) -Z<(X,z)Z>(z,X)] . (44) 



The leading contribution to the self energies Y1h u (x, y) and TjH d {x,y) from the inter- 
actions of Eq. ff22l are (see Figure [3]) 

^H u (x,y) = g(x,y)G HL (x,y) + h(x,y)G H > d (x,y), 

t Hd (x,y) = g(x,y)G H '(x,y) + h(x,y)G H i(x,y), (45) 



14 



where 



g(x,y) = m 2 uu ,{x)m 2 * uu ,{y) = [c* 2 ^T){x) + ci/i'*E(a:)] [c 2 //*S(y) + c^'E^y)] 

= m^(x)m 2 ^(j,) = [A t E(x) + ^(x)] [A{tf (y) + A'^y)] 

g(x,y) = m 2 dd ,{x)m 2 * dd ,{y) = [c 2 ^'*S(x) + c*fiE\x)] [(^(i'E\y) + cifi*E(y)] 

h(x,y) = m 2 u ,(x)m 2 * du ,(y) = [A 2 Z(x) + A' 2 ^ (x)] [A* 2 tf(y) + A%E(y)] (46) 

The CP- violating part of Eq. (|44|) . evaluated for H u , is the imaginary part of Green 
functions, 



1h u 



Xo 



i dz dz < [g(x, z) - g(z, x)] Im G^(x,z)G^ u (z,x) - G^(x,z)G^ u (z,x) 



+ [h(x,z) - h(z, x)] Im Gj» (x,z)G^ u (z,x) - G|, (x, z)G^ u (z,x) 



(47) 



We can express the scalar Green's function in terms of the Bose-Einstein distribution and 
the spectral density of a scalar field 



G*{x,y) 



d A k 



where the equilibrium distribution functions are: 



(48) 



(49) 
(50) 



with ub(x) = l/[exp(x/T) — 1]. The spectral density in the limit of small decay width is 



p(k ,k) 



1 



1 



(k° + ie + iV) 2 - u 2 (\k\) (jfc° -is- iT) 2 - uo 2 {\k\) 



(51) 



where u; 2 (|/c|) = k 2 + M 2 . The mass M and width V are the thermal mass and width, 
respectively. 

The coefficients [g(x, z) — g(z, x)] and [h(x, z) — h(z, x)] can be calculated to the first 
order in the expansion about x — z. 



9{x,z) - g{z,x) 



2% [|ci/i'| 2 — |c2/x| 2 ] sin [9(x) — 9(z)} u(x)u(z) 
2i [Id/if - \c 2 fi\ 2 } (x - zY [d^B{x)\ u 2 {x) + ■ 



(52) 



When these are inserted in Eq. (T4"T|) . the fact that the spectral density is isotropic in space, 
implies that only the time component is non- vanishing, and we can make the replacements, 

g(x,z) -g(z,x) = 2i [|ci/i| 2 - |c 2y u| 2 ] (x - z)°] ( ^-v w J u 2 (x), (53) 



h(x, z) — h(z, x) 



u 2 (x) 



(54) 
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and thus, 
where 

T H U H{ 



(fS) U2{X) { ^ " |C2/i|2 ] 1h -< + ^ ~ ^ X ^ H '} ' (55) 



dk- 



k 2 



x 



1 + 2Re[n B (c^ H ; + ir H ()])I(w Hu , r Hu , w H ;, T R > 



o 2ti 2 uj h ,uj Hu 
+ (1 + 2Re[n B (^ Hu + iT H J])I(^, r H (, w Hu , r H J 

-2(Im[n B (wH u + lThu)] + Im[n B (^ H ; + ir H ;)])G(u;H u , T HuJ uj r >, r H ;) } , (56) 



and the functions / and G are given by 



I(a, b, c, d) 



1 



1 



2[(a + c) 2 + (6 + d) 2 " 
1 



sin 



2[(a-c) 2 + (6 + d) 2 ] 



sin 



2 arctan 
2 arctan 



a + c 
b + d 
a — c 
b + d 



(57) 



G(a, b, c, d) 



1 



1 



2 [(a + c) 2 + (6 + d) 2 } 
1 1 



cos 



2[(a-c) 2 + (6 + rf) 2 ] 



cos 



2 arctan 
2 arctan 



a + c 
b + d 
a — c 



b + d 



In exactly the same way, we derive 



(5f 



: (l^) M2(X) { ^ l/i|2 " |C2//|2 ] J ^ + [ |Aa|2 " ]A ' 2 ^ Xh ^} ■ (59) 

The over-all magnitude of the CP- violating source depends sensitively on several param- 
eters which have up until now not played a large role in deriving our ressults. Thus, we 
content ourselves with the order of magnitude estimate based on the sample parameters 
for the a potential and bubble wall velocity and profile, and assuming the thermal masses 
and widths are roughly TeV, and that the H-H-H 1 dimensionful interaction^ are of order 
TeV. Evaluating the integrals numerically, we find 7^ ~ 0.01 TeV A6/L w v w u 2 , which for 
the choice of sample parameters described above leads to 7^ ~ 10 9 GeV 4 , and acquires 
non-vanishing values only at the position of the bubble wall, where the Higgs fields are 
varying. 



5 The thermal widthes of the Higgs and spectator Higgs are dominated by the trilinear interactions 
E-H-H'. 
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The results presented above rely on a method similar to the one used in Ref. [19] 
in the Minimal Supersymmetric Standard Model case. Further improvements to this 
method were performed in Refs. [21] . by considering the all-order resummation of the 
Higgs mass insertion efffects. The result was a mild suppression of the results in the case 
of degenerate masses, as well as new relevant terms away from the degenerate mass regime. 
Furthermore, in Ref. [22], based on self-consistency arguments, a more detailed analysis 
of the relation between the CP-violating sources and the currents induced by the Higgs 
fields was performed, leading to the presence of higher order derivatives in the sources, as 
first suggested in Ref. [17]. The final result of these investigations is a suppression of the 
baryon asymmetry by a factor of a few in the degenerate mass regime compared to the 
one obtained in Ref. [19], as well as a power- law suppression with the masses when they 
move away from the degenerate region. 

One of the weaknesses of the above described work is the lack of a rigorous derivation 
of the sources for the diffusion equation. A derivation of the semiclassical forces in the 
transport equations by means of the dynamics of the two-point function in the Schwinger- 
Keldysh formalism was performed in Refs. [231 [2H [221 [2S] • In particular, in Ref. [2S] a 
consistent treatment of the fermion mixing effects was performed, and a careful derivation 
of the sources in the diffusion equations was obtained. The final result was smaller 
by a factor of a few to an order of magnitude than the result obtained in Ref. [22] in 
the degenerate mass regime, and a stronger than power-law suppression away from the 
degenerate case. In our work we are interested in an order of magnitude estimate of 
the result for the baryon asymmetry, and for that purpose the results of this section are 
sufficient. However, while an enhancement of the sources by a factor of a few would be 
easy to obtain by a careful choice of the parameters of our model, an enhancement by 
several order of magnitudes woulde prove very difficult. For these reasons, it would be 
interesting to pursue a more general treatment using the techniques of Ref. [26J to make 
a detailed exploration of the full parameter space consistent with the observed baryon 
asymmetry. 



3.5 Results 

We assemble the results and make a prediction for the density of lepton number stored 
in the third family, L3. We have solved the diffusion equations in two ways; the first is 
a brute force numerical solution of the coupled differential equations, while the second 
proceeds by making some simplifying approximations which allow us to determine an 
analytic solution to the diffusion equations. Our analytic solution simplifies the problem 
by assuming that T y , Tq C d and (for z < 0) r x are all strong enough that they enforce 
near-equilibrium relations among the number densities of the species participating in the 
interactions. Thus, we have, 

Q h t 1 

0, 



k-Q kh k t T y 
Q t b 1 
kQ k t k b Iqcd 
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and for z < 0, 



3 ^ + ^l 
k Q k L T 1 



0. 



(61) 



Following the treatment of the usual EW case |37J, we take linear combinations of the dif- 
fusion equations (1391) which are independent of T y , Tqcd, and I\ and use the equilibrium 
relations to express the remaining two equations in terms of Q and h. The result can be 
expressed matrix differential equation, 



M 



' Q" ' 


+ N 


Q' 








h" 




h! 








(62) 



where we simplfy 7^ as constant for — L w < z < and M and N are constant matrices, 



functions of the .D's, fc's, and v w , 
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= v w 






41 


1 




9 
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(63) 



We convert Eq. f)62p into a first order differential equation by defining \1/ = [Q , h] T , 
= We solve the equations seperately for z < L w (region I), — L w < z < (region 
II), and z > (region III) and match ty and H across the boundaries. 

In regions I and III, 7^ = 0, and the solutions are those of the corresponding homo- 
geneous equation, 



H lj m(z) = H H (z)V l , m = Exp[-M^Nz]V ljm 

tyi,m(z) = / dx H H (x) Pyn + Ci, 

J —00 



(64) 
(65) 



where Cini and Pim are vectors specifying the boundary conditions. The exponential 
terms grow with z, and thus the requirement that ty(z — > —00) — > requires C\ = and 
the requirement that ty remain finite as z — > +00 requires Z>m = 0. 

In region II, the solutions will be the sum of a homogeneous solution with integration 
constants C\\ and T>u, and a particular solution, 



So, 



H11 (^) 
*n (z) 



-H H {z) 



in 




H„(z) V n + Hp(z) 



dx {H H (x) V u + Hp (x)} + C n 



(66) 

(67) 
(68) 
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and matching this to the solutions in regions I and II determines T>\ = T>n, and, 



o 



dxHjj 1 (x) M~ l 



dxHn (x) T>n 



lh 




C 



in 



;-0 px 

dxH H (x) V n - dx dyH H (x) H H X (y) M" 1 



lh 





(69) 
(70) 
(71) 



Note that ^(z) = Cm for z > 0, so this last expression is in fact the final densities 
produced by the phase transition. 

Since L w is small, we can derive an approximate form based on the limit L w — > 0. To 
leading order in L w , we have 



C 



in 



dxH H (x) Vu = -L w N~ l 



lh 




(72) 



where 



V 



ii 



L„,M _1 



lh 




(73) 



Using the expression for iV 1 , we find the final lepton number density 



-Q 



" 227v v , 



(74) 



Recalling equation (j581) that ^ h oc v w /L w , we find that our final result is approximately 
independent of the diffusion constants, the bubble wall velocity, and the bubble wall 
width, as long as Ti is fast enough. 

Assembling the results, in Figure H] we plot the densities normalized to the entropy, 
s ~ 27r 2 /45g*T 3 , where ~ 100. The densities L 3 , t, and b are determined using the 
relations Eqs. (I60H61I) . The wall region is too small (2 x 10~ 3 GeV -1 ), and is only in the 
figure For the chosen parameters, the ratio of the lepton density to the entropy density 
at is given by 

L 3 /s ~ lO" 4 , (75) 

which results, after the inclusion of the dilution factor, Eq. f[T2T) in a final baryon asym- 
metry after the electroweak phase transition of approximately 



^ ^ io- 10 , 



(76) 
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Figure 4: Particle number densities normalized to entropy as a function of spatial position 
z for a bubble whose wall is at z ~ and parameters as described in the text. From top 
to bottom, the curves represent the following densities: h, Q3, bn, L 3 , and £r. 
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exactly as observed. Of course, the particular value is highly dependent on our choice of 
parameters, but the ability of the Topflavor model to produce this value is not; the fact 
that the order of magnitude comes out correctly is indicitive of the fact that for natural 
values of parameters, this model can produce an appropriate baryon number. 

From our numerical integration of the differential equations, we can relax our equilib- 
rium assumptions and examine the effect of finite values for the rates, particularly r\. In 
Figure Owe present the final lepton density as a function of Ti (artificially assuming that 
the CP- violating source, bubble parameters, and critical temperature are unchanged). 
For small Fx, we see that the final lepton asymmetry is linearly proportional to Fx- In 
this regime, the generated Higgs/top number densities are a constant background source, 
only a small fraction of which is converted into leptons by Fx- For large Ti instead, 
the equilibrium condition is reached and the dependence on Fx saturates. We see from 
Figure [5] that our chosen parameters are just at the turn-on of this saturation region. 



4 Discussion 

We have explored the possibility that the baryon asymmetry of the Universe might have 
come from new gauge dynamics at the TeV scale. The primary impediment to realizing 
this idea is the washout by ordinary EW sphalerons, which we have avoided by having new 
physics which coupled asymmetrically to the three (lepton) generations^. In particular, 
we showed that the Topflavor model, with SU(2)x x S77(2) 2 breaking at the TeV scale, 
can easily produce an appropriate baryon asymmetry for natural values of its parameters. 

It is interesting to compare the results obtained in Eq. (1751) with the much smaller 
ones that are obtained by similar methods in the MSSM [19J. The difference stems from 
three main different factors: First, the magnitude of Fx is much larger than the EW 
sphaleron rate r ws ~ 10~ 4 GeV. As shown in Fig. [51 this induces an enhancement of a 
few orders of magnitude in the final result. Second, the value of A6 ~ 0.2 is two orders 
of magnitude larger than the typical values of A/3 ~ O (10 -3 ) obtained in the MSSM for 
moderate values of the CP-odd Higgs mass. Third, in the MSSM there is a suppression 
proceeding from the effective cancellation of chiral charges discussed in Refs. [5U |37] that 
is not present in the model analyzed in this article. 

We expect that the LHC can extensively test this idea. First, by discovering the 
SM Higgs boson and searching for the existence of other elements strongly coupled to it, 
the LHC can help construct the picture of the ordinary EW phase transition as first or 
second order. This would, for example, confirm or rule out the MSSM scenario of EW 
baryogenesis. Second, by discovering new elements, in the case of Topflavor, the heavy W 
and Z' bosons with masses of about 2 TeV and preferred decays into third family fermions, 
should be easily discovered at the LHC [55]. Finally, by determining the properties of 
the new elements, one may determine whether or not they are viable as a mechanism of 

6 A parallel idea, which we chose not to explore in detail, is to have new physics coupled asymmetrically 
between quarks and leptons, as in [53] . Anomaly cancellation is very challenging in this framework, and 
the resulting models are therefore somewhat contrived. 
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r 1 (GeV) 

Figure 5: Lepton number density normalized to entropy as a function of F±. The red 
dashed line is the bound on Ti from non-observation of proton decay mediated by instan- 
tons, whereas the blue dotted line is the bound inferred by requiring that the broken phase 
minimum remain the true vacuum at T = 0. The star indicates the sample parameter 
point considered in the text. 
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baryogenesis. In the case of Topflavor, the need for a first order phase transition in the 
breakdown of SU(2) 1 x SU (2) 2 is connected to both a small quartic for the Higgs a and a 
reasonably strong gauge coupling g\. The first of these has the dramatic consequence that 
the mass of a is well below the symmetry-breaking scale, m 2 a ~ Am 2 . For the parameters 
we explore in detail, m a ~ 600 GeV, making a the lightest of the new particles. It couples 
strongly to the Higgs sector, and to the new gauge bosons. 

An extended gauge sector at the TeV scale has motivations from many mysteries in 
particle physics. In addition, it allows a much less constrained structure to explore the 
idea of baryogenesis through gauge symmetry breaking. The SM picture arising from the 
EW breakdown is economical, but is under assault from the null LEP Higgs searches, 
favoring a second order phase transition, and from the lack of sufficient CP violation to 
produce enough asymmetry. Alternatives based on the ordinary EW transition such as in 
the MSSM can be viable, but remain constrained by the Higgs mass, and from the fact 
that the new sources of CP violation remain tightly connected to low energy CP- violating 
observables. Extended gauge symmetries can naturally have first order transitions, and 
CP violation is further removed from low energy observables. The primary obstacle is 
the fact that the EW sphalerons will try to erase any generated baryon asymmetry with 
B — L = 0, but this is overcome if the SM matter does not couple universally to the new 
dynamics. Ultimately, the LHC will help in resolving this question. 



5 ACKNOWLEDGMENTS 

We would like to thank Junhai Kang, and David Morrisey for helpful conversations and 
Ignacy Sawicki for computing assistance. Work at ANL is supported in part by the US 
DOE, Div. of HEP, Contract W-31-109-ENG-38. J.S. would like to acknowledge Prospects 
in Theoretical Physics 05 for part of the initial work was done there, T.T. and C.W. would 
like to thank the Aspen Center for Physics were an important part of this work has been 
performed. 



References 

[1] A. D. Sakharov, Pis'ma Zh. Eksp. Teor. Fiz. 5, 32 (1967). 

[2] M. Yoshimura, Phys. Rev. Lett. 41, 281 (1978) [Erratum-ibid. 42, 746 (1979)]. 

[3] S. Dimopoulos and L. Susskind, Phys. Rev. D 18, 4500 (1978). 

[4] S. Weinberg, Phys. Rev. Lett. 42, 850 (1979). 

[5] E. W. Kolb and M. S. Turner, The Early Universe, Addison- Wesley Publishing Com- 
pany (1990). 

[6] M. Fukugita and T. Yanagida, Phys. Lett. B 174, 45 (1986). 



23 



[7] W. Buchmuller and M. Plumacher, Int. J. Mod. Phys. A 15, 5047 (2000), and refer- 
ences therein. 



[9 
[10 

[11 

[12 
[13 
[14 
[15 
[16 
[17 
[18 
[19 

[20 
[21 

[22 
[23 

[24 

[25 

[26 



I. Affleck and M. Dine, Nucl. Phys. B 249, 361 (1985). 

V. A. Kuzmin, V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B 155, 36 (1985). 

For reviews, see, M. Trodden, Rev. Mod. Phys. 71, 1463 (1999); V.A. Rubakov, M. 
E. Shaposhnikov, Usp. Fiz.Nauk 166 493 (1996); Phys. Usp. 39 461 (1996); A. G. 
Cohen, D. B. Kaplan, A. E. Nelson, Ann. Rev. Nucl. Part. Sci. 43 27 (1993). 

R. Barate et al. [LEP Working Group for Higgs boson searches], Phys. Lett. B 565, 
61 (2003). 

M. Pietroni, Nucl. Phys. B 402, 27 (1993). 

A. T. Davies, C. D. Froggatt and R. G. Moorhouse, Phys. Lett. B 372, 88 (1996). 

S. J. Huber and M. G. Schmidt, Nucl. Phys. B 606, 183 (2001). 

J. Kang, P. Langacker, T. j. Li and T. Liu, Phys. Rev. Lett. 94, 061801 (2005). 

A. Menon, D. E. Morrissey and C. E. M. Wagner, Phys. Rev. D 70, 035005 (2004). 

J. M. Cline, M. Joyce and K. Kainulainen, JHEP 0007, 018 (2000). 

J. M. Cline and K. Kainulainen, Phys. Rev. Lett. 85, 5519 (2000). 

M. Carena, M. Quiros, A. Riotto, I. Vilja and C. E. M. Wagner, Nucl. Phys. B 503, 
387 (1997). 

M. Carena, M. Quiros and C. E. M. Wagner, Nucl. Phys. B 524, 3 (1998). 

M. Carena, J. M. Moreno, M. Quiros, M. Seco and C. E. M. Wagner, Nucl. Phys. B 
599, 158 (2001). 

M. Carena, M. Quiros, M. Seco and C. E. M. Wagner, Nucl. Phys. B 650, 24 (2003). 

K. Kainulainen, T. Prokopec, M. G. Schmidt and S. Weinstock, JHEP 0106, 031 
(2001). 

K. Kainulainen, T. Prokopec, M. G. Schmidt and S. Weinstock, Phys. Rev. D 66, 
043502 (2002). 

T. Prokopec, M. G. S chmidt and S. Weinstock, Annals Phys. 314, 208 (2004) 
|arXiv:hep-ph/0312110|; T. Prokopec, M. G. Schmidt and S. Weinstock, Annals Phys. 
314, 267 (2004) [arXiv:hep-ph/0406i40| . 

T. Ko nstandin, T. Prokopec M. G. Schmidt and M. Seco, Nucl. Phys. B 738, 1 
(2006) [arXiv:hep-ph/0505103| . 

24 



[27; 

[28 
[29 
[30 
[31 

[32 
[33 

[34; 

[35 



[36 

[37; 

[38 
[39 

[40 
[41 
[42 
[43 
[44; 
[45 
[46 
[47 



M. Joyce, Phys. Rev. D 55, 1875 (1997). T. Prokopec, Phys. Lett. B 483, 1 (2000). 

M. Joyce and T. Prokopec, Phys. Rev. D 57, 6022 (1998). 

S. Davidson, M. Losada and A. Riotto, Phys. Rev. Lett. 84, 4284 (2000). 

G. Servant, JHEP 0201, 044 (2002). 

M. Carena, A. Megevand, M. Quiros and C. E. M. Wagner, Nucl. Phys. B 716, 319 
(2005). 

A. Provenza, M. Quiros and P. Ullio, JHEP 0510, 048 (2005). 

V. A. Kuzmin, V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B 191, 171 
(1987). 

H. K. Dreiner and G. G. Ross, Nucl. Phys. B 410, 188 (1993). 

R. S. Chivukula, E. H. Simmons and J. Terning, Phys. Rev. D 53, 5258 (1996); 
E. Malkawi, T. Tait and C. P. Yuan, Phys. Lett. B 385, 304 (1996); D. J. Muller 
and S. Nandi, Phys. Lett. B 383, 345 (1996); P. Batra, A. Delgado, D. E. Kaplan 
and T. M. P. Tait, JHEP 0402, 043 (2004). 

D. E. Morrissey, T. M. P. Tait and C. E. M. Wagner, Phys. Rev. D 72, 095003 (2005). 

P. Huet and A. E. Nelson, Phys. Rev. D 53, 4578 (1996). 

S. Coleman, V. Glaser and A. Martin, Commun. Math. Phys. 58 211 (1978). 

G. D. Moore and T. Prokopec, Phys. Rev. D 52, 7182 (1995) |arXiv:hep-ph/9506475| ; 
G. D. Moore and T. Prokopec, Phys. Rev. Lett. 75, 777 (1995) 
|arXiv:hep-ph/9503296l . 

A. Kosowsky, M. S. Turner and R. Watkins, Phys. Rev. Lett. 69, 2026 (1992). 

R. Apreda, M. Maggiore, A. Nicolis and A. Riotto, Nucl. Phys. B 631, 342 (2002). 

A. Nicolis, Class. Quant. Grav. 21, L27 (2004). 

C. Grojean and G. Servant, |arXiv:hep-ph/0607107| . 

L. Randall and G. Servant, |arXiv:hep-ph/0607158] . 

J. Shu, T. M. P. Tait and C. E. M. Wagner, in progress. 

A. D. Linde, Nucl. Phys. B 216, 421 (1983) [Erratum-ibid. B 223, 544 (1983)]. 

M. Dine, R. G. Leigh, P. Y. Huet, A. D. Linde and D. A. Linde, Phys. Rev. D 46, 
550 (1992). 



[48] A. G. Cohen, D. B. Kaplan and A. E. Nelson, Phys. Lett. B 336, 41 (1994). 

25 



[49] M. Joyce, T. Prokopec and N. Turok, Phys. Rev. D 53, 2930 (1996). 

[50] P. Arnold, D. Son and L. G. Yaffe, Phys. Rev. D 55, 6264 (1997); D. Bodeker, 
G. D. Moore and K. Rummukainen , Nucl. Phys. Proc. Suppl. 83, 583 (2000); 
G. D. Moore, larXiv:hep-ph/0009161| . 

[51] A. Riotto, Phys. Rev. D 53, 5834 (1996). 

[52] A. Riotto, Nucl. Phys. B 518, 339 (1998). 

[53] H. Georgi, E. Jenkins and E. H. Simmons, Nucl. Phys. B 331, 541 (1990). 
[54] G. F. Giudice and M. E. Shaposhnikov, Phys. Lett. B 326, 118 (1994). 



[55] Z. Sullivan, arXiv:hep-ph/0306266|. 



26 



